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Abstract 



In this paper dynamic von Karman equations with localized interior damping supported in a 
boundary collar are considered. Hadamard well-posedness for von Karman plates with various 
types of nonlinear damping are well-known, and the long-time behavior of nonlinear plates has 
been a topic of recent interest. Since the von Karman plate system is of "hyperbolic type" with 
critical nonlinearity (noncompact with respect to the phase space), this latter topic is particularly 
challenging in the case of geometrically constrained and nonlinear damping. In this paper we first 
show the existence of a compact global attractor for finite-energy solutions, and we then prove that 
the attractor is both smooth and finite dimensional. Thus, the hyperbolic-like flow is stabilized 
asymptotically to a smooth and finite dimensional set. 

Key terms: dynamical systems, long-time behavior, global attractors, nonlinear plates, nonlinear 
damping, localized damping 

1 Introduction 

We consider the evolution of a nonlinear von Karman plate subject to nonlinear frictional damping 
with essential support in a boundary collar. Our aim is to consider the long-time behavior of the 
corresponding evolution. This includes studying (a) existence of global attractor which captures long- 
time behavior of the dynamics, and (b) properties of this attractor, such as smoothness and finite 
dimensionality. 

In short, our goal is to show that the original infinite dimensional and non-smooth dynamics of 
hyperbolic type can be reduced (asymptotically) to a finite dimensional and regular set, with respect 
to the topology of "finite energy" . The latter is associated with weak (or generalized) solutions of the 
underlying semigroup for the dynamics. This type of result then allows the implementation of tools 
from finite dimensional control theory in order to achieve a preassigned outcome for the dynamics. 
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1.1 Model and Energies 

Let ft C M 2 be a bounded domain with dft = T taken to be sufficiently smooth. We consider a plate 
model where the real-valued function u(x,y;t) models the out-of-plane displacement of a plate with 
negligible thickness. Then the von Karman model [T^lIM] requires that u satisfies 

u tt + A 2 m + d(x)g(ut) = fv{u) + p in ft x (0, oo) = Q (1.1) 
u\ t=Q = u , u t \ t ^ Q = Ul . 

The von Karman nonlinearity 

f v =[v(u)+F ,u] (1.2) 
is given in terms of (a) the Airy Stress function v(u), satisfying 

A 2 v(u) = ~[u,u] in fl (1.3) 
d v v(u) — v{u) = on r, 

and (b) the von Karman bracket given by 

[U, W] = U XX Wyy + UyyW XX ~ 2U xy W X y. (l.4) 

The internal force Fq £ H (fl) n Hq(£1), 9 > 3, and external force p £ £2^) play a n essential 
role in shaping the nontrivial stationary solutions. (In this paper H 8 (D) denotes the Sobolev space of 
order seRon domain D.) In the absence of these forces, the stationary solution of the corresponding 
nonlinear boundary value problem becomes trivial and simply reduces to zero. 

In this treatment we focus on the stabilizing properties of the damping term d{n)g{u t ). In particu- 
lar, we take <?(•) £ C(R) to be a monotone increasing function, with g(0) = and further boundedness 
and smoothness assumptions to be imposed later; additionally, d(x) = d w (x) is a nonnegative i^fi) 
localizing function which restricts the damping term g(ut) to a particular subset w C fi. This is to 
say u> C supp d or d(x) > Co > for x £ u>. Initially we will take us to be a general set oj CC 57, but 
more specifically, we are interested in taking uj to be an open collar of the boundary T. This type of 
damping represents localized, viscous damping near the boundary T. 

The boundary conditions we consider for the plate are: 

1. Clamped, denoted (C) 

u = d v u = in rx(0,oo) = E. (1.5) 

2. Hinged (simply-supported), which we denote by (H) 

u = Au = in S. (1.6) 
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3. Free-type, denoted by (F) 



B x u = Au + (1 - (i)B x u = on T u 
B 2 u = d v Au + (1 — \i)B 2 u — — fiv? — on Ti, 

u = d v u = (clamped) on r , (1.7) 

where we have partitioned the boundary r = To U Ti (with To possibly empty). For simplicity 
we assume that To H Ti = 0. Otherwise the regularity theory for elliptic problems with mixed 
boundary conditions must be invoked. The boundary operators B\ and B 2 are given by |34j : 

Biu = 2viv 2 u xy - v\u yy - v 2 u xx , 

B 2 u = d T [{yl - u xy + v x v 2 (u vy - u xx )] , 

where v = (i>x, v 2 ) is the outer normal to T, r — {—v 2 , v\) is the unit tangent vector along T. The 
parameters fix and /3 are nonnegative, the constant < fj, < 1 has the meaning of the Poisson 
modulus. 



Notation: Note, when referencing the plate equation above in (1.1 1 wc will write ([Ll])(C), ([l~T|)(H) 



or (1.1 )(F) to indicate which boundary conditions we are taking. We write the norm in H S {D) as 
|| • || s and || • 1 1 = II ' ||l 2 (£>)) f° r simplicity (when the meaning is clear from context) norms and inner 
products written without subscript ( (•,•), || • || ), are taken to be L 2 (D) of the appropriate domain D. 
Additionally, we employ the notation that Hq(D) gives the closure of Cg°(D) in the || • || s norm. 

The von Karman plate equation is well-known in nonlinear elasticity, and constitutes a basic model 
to describe the nonlinear oscillations of a thin plate with large displacements |34j (and references 
therein). In particular, we take the thickness of the plate to be negligible (as is usual in the modeling 
of thin structures [H]). 

Remark 1.1. It is worth noting that the von Karman plate model can accomodate plates with non- 
negligible thickness - the equation in ( |1.1[ ) then gives the vertical displacement of the central plane 
of the plate. This is tantamount to adding the term —^Au tt , 7 > to the LHS of (JOJ. This term 
corresponds to rotational inertia in the filaments of the plate, and (a) is regularizing from the energetic 
point of view and (b) forces the dynamics of the plate to be hyperbolic. In this treatment we take 
7 = 0, since it constitutes the most challenging problem mathematically, however, a future manuscript 
will address the case 7 > and the limiting problem (convergence of solutions and attractors) as 
7\0. 

The energies associated to the above equation are given by (in the case of clamped (C) or hinged 
(H) boundary conditions) 

E(t)=^(\\Au\f + \\u t \\ 2 ), 

E{t) =E{t) + \\\Av{uW 
#(t) =E(t)+U(u), 
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where 



Il(u) = -I (\Av(u)\ 2 -2[F ,u}u-4pu). (1.8) 



v. 



The above (linear) energy E(t) dictates our state space T-L, which depends on boundary conditions. 
In the case of clamped boundary conditions (C) we have Hi = H 2 (£l) x £2(0). For hinged boundary 
conditions (H) we have U 2 = (H 2 n Ho)(0) x L 2 (Cl). 

Lastly, for free boundary conditions (F) we have H3 = (H 2 n H^ Ta )(Cl) x L 2 (fl) (where H 2 r (fi) 
is the Sobolev space H 2 (VL) with clamped conditions on To); the potential energy in this case is given 
by the bilinear form 



a(u,v) — / a (it, v) + fii / mw, (1-9) 
Jn JTi 

where 

ffl(U, V) = U XX V XX + UyyVyy + [l(U XX Vyy + UyyV XX ) + 2(l - [l)U X yV X y. (^lO) 



Then the energy becomes 



The total energy becomes 



%) = £( t ) + I||A<;(«)|| 2 + f^ U 4 dr 



^(t) = E(t) + n(«(t)) + i)9^ u *( t ). 



Remark 1.2. We note that this last form of the energy described by the bilinear form a(u, v) can also 
be applied to clamped or hinged boundary conditions. Indeed, in this latter case the bilinear form 
a(u,u) collapses just to ||Au|| 2 . 

It will be convenient to introduce an elliptic operator A: @(A) C L 2 (fl) — > L 2 (£l) given by Au = 
A 2 u, where 2>(A) incorporates the corresponding boundary conditions (clamped, hinged, or free). It 
is useful to note that by elliptic regularity 

{H$(Sl) clamped BC 

{H 2 nH^)(n) hinged BC 
(H 2 H H% iT0 )(O) freeBC 

It is important to note the total potential energy may not be positive, or even not bounded from 
below. This is due to the presence of internal force Fq which may drive the energy to —00. However, 
the presence of the von Karman bracket in the model, along with appropriate regularity properties 
imposed on Fq, assures that the energy is bounded from below. This can be seen from the following 
lemma [HQS: 



Lemma 1.1. Let u € ^(A 1 / 2 ), p e L 2 (Q), and F e H^(Q) n H e {Vt), 9 > 3. Then, V e > th 
exists M(e, \\p\\, \\FoW0) = M c>Pi f < 00 such that in the clamped and hinged case 



ere 
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<e(||^/ 2 n|| 2 + ||Ai;(u)|| 2 )+M e)PjFo 
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and in the free case with /3 > 0, 



IMI 2 < eiWA^uW 2 + \\Av(uW + f IHIi 4( r)) + M *M 

As a consequence we have the following bounds from below for the energy: 
There exist positive constants to, c, M, C such that 

-m + cE(t) <<?(*) < M + CE(t) (1.11) 
-m + cE(t)<£(t)<h(E(t)) (1.12) 

where /i(s) denotes a continuous function. 
1.2 Motivation and Literature 

Well-posedness for von Karman's plate equation with interior and/or boundary dissipation has been 
known for some time for smooth solutions in the case of homogeneous [5] or inhomogeneous nonlinear 
boundary conditions (TSJ, |2T] and references therein. The issue of well-posedness for 'weak' (finite- 
energy) solutions is more recent |12U21| . In this paper, we are interested in homogeneous type boundary 
conditions and we will be considering generalized nonlinear semigroup solutions |H [33] which also can 
be shown to be weak variational solutions. For a detailed and complete discussion regarding the 
wellposedness and regularity of von Karman solutions the reader is referred to [12 [25] . In the context 
of this paper we will need the following well-posedness result, which is contingent upon the recently 



shown sharp regularity of the Airy Stress function in (1.3) [2"Tl 112) : 



Theorem 1.2. With reference to problem \l.l\ C) with initial data (uq,u\) £ Hi, or \l.l\ H) with initial 
data (uq,ui) £ H2, or \l.l\ F) with intial data (u$,ui) £ H3, there exists a unique global solution of 
finite-energy (i.e. (u,u t ) £ C([0,T];'Hi) fori = 1, 2, 3 resp., for any T > 0). Additionally, (u,Ut) 
depends continuously on (uq,ui) £ Hi. 

Thus, for any initial data in the finite energy space (u , ui) € H., there exists a well defined scmiflow 
(nonlinear semigroup) St{u^,ui) = (u(t),u t (t)^ £ H which varies continuously with respect to the 
initial data in T~L. The domain of the corresponding generator A(u, v) = (v,—Au—d(yL)g{v)-\-fv{u)+p^j 

is given by 9{A) = {(u,v) € 2>(A 1/2 ) x @(A 1/2 );Au + d(x)g(v) £ L 2 (f2)}. For initial data taken in 
0(A), the corresponding solutions are regular and remain invariant in 0(A) [H EM ES]- With an 
additional assumption that g(s) is bounded polynomially at infinity, one has 0(A) C H 4 (fl) x H 2 (il). 
Equipped with the regularity of the domain 0(A), one derives the energy identity for all regular 
solutions. Due to the density of the embedding 0(A) C %, monotonicity of the damping, and sharp 



regularity of the Airy stress function (see Lemma 3.4 1 the same energy equality remains valid for all 



generalized solutions corresponding to any boundary conditions under consideration. Thus we have 
the energy identity for boundary conditions (C), (H), or (F) satisfied for all generalized (semigroup) 
solutions (complete details of this argument are given in |12|). 



This equality reads: for all < s < t, strong and generalized solutions u to (1.1) satisfy 



S(t) + / / d(x)g(u t )u t = S(s) (1.13) 
n 



5 



With the well-posedness of the semiflow established in Theorem 1.1, it is natural to investigate long 



time behavior of the dynamical system generated by (1.1). It is clear from (1.131 that the essential 
mechanism for dissipating the energy is the damping term d{x)g(u t ). In the simplest possible scenario 
when p — Fq = the energy function $(t) is equivalent topologically to the norm of the phase 
space %. Since <?(t) is nonincreasing on the trajectories, it becomes a Lyapunov function for the 
corresponding nonlinear dynamical system, whose only equilibrium is the zero point. If one assumes 
that d(x) > 0, a.e.in Q, it is well known that S'it) becomes a strict Lyapunov function and zero 
equilibrium is strongly stable. However, the above condition imposed on d{x) is not sufficient to 
guarantee uniform convergence to the equilibrium (this is also the case for linear dynamics without 
the von Karman term). In order to secure uniform convergence or, more generally, convergence to 
a compact attractor, a stronger form of the damping is necessary. For example, d(x) > cq > 0, 
x G £1 and g(s) = as, a > 0, provides a classical model for which uniform convergence to zero in the 
absence of external/internal forcing (or more generally to an attractor) can be shown HI [T2J. 134] . 
The goal in this paper is to consider nonlinear damping of a reduced essential support whereby the 
inequality d(x) > Co > will be enforced only in a small set w CC !J, while the dynamics will be 
forced by nontrivial sources p, Fq. Existence of a compact and possibly smooth finite dimensional 
attracting set for the dynamics generated by ( |1.1[ ) with boundary conditions (C), (H), or (F) and 
geometrically constrained dissipation is of great physical interest. Such a result is tantamount to 
asserting that the infinite dimensional, non-smooth dynamics are asymptotically reduced to a smooth 
and finite dimensional set. While such a reduction is expected for dynamical systems that exhibit some 
smoothing effects (e.g. parabolic-like) [HI (37J El HS 03 GM US] , it is a much less evident phenomena 
in the case of hyperbolic-like dynamics, where the 'taking-off ' of the dynamics produces no smoothing 
effect. The role of the frictional damping in such a system is instrumental; in fact, it is the induced 
friction that creates a stabilizing and asymptotically regularizing effect on the evolution, ultimately 
reducing it to a compact set. On the other hand it is well known that the hyperbolic-like dynamics 
can not be stabilized by a compact feedback operator [30] (and references therein). This is due to 
the fact that instabilities in the system are inherently infinite dimensional and the essential part of 
the spectrum can not be dislodged by a compact perturbation. Thus, any effective damping cannot 
be compact (with respect to the phase space). The above feature combined with (a) nonlinearity of 
the damping and (b) lack of compactness of the nonlinear von Karman source makes the analysis 
of long-time behavior for this class of systems challenging. In fact, critical exponent nonlinearities 
and nonlinear dissipation are known to constitute endemic difficulties in the study of hyperbolic-like 
systems [22] . 

To orient the reader and to provide some perspective for the problem studied, we shall briefly 
describe some of the principal contributions to this area of research. A detailed account is given in 

In the discussion of global attractors for von Karman evolution equations, we must distinguish 
between two types of dynamics for the problem: (a) the rotational case (as addressed above) when 



the term —~/Au t t, 7 > is added to the LHS of (1.1 1 and (b) nonrotational (7 = 0). In case (a), we 



note that the von Karman nonlinearity (in the finite energy topology) is compact, which considerably 
simplifies the analysis of long-time dynamics. In the latter case (b) (which we consider here), a very 
different type of analysis is needed. Here, we shall focus on part (b) only. In fact, the very first 
contribution to this problem is a pioneering paper [5] where the existence of weak attractors with a 
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linear, fully supported damping was demonstrated. Later on, owing to new results on the regularity of 
Airy's stress function [211 H2], weak attractors were improved to strong attractors, and the restriction 
of linear damping was removed in order to allow nonlinear, monotone damping [11] . In order to 
incorporate fully nonlinear interior damping, [11] assumes that the dissipation parameter is sufficiently 
large. This restriction was later removed in [27] , whose paper introduces a very clever way of bypassing 
a lack of compactness and replacing it with an "iterated convergence" trick. Further studies of the 
attractor (including properties such as dimensionality and smoothness) in the fully nonlinear setup, 
without "size" restrictions imposed on the parameters, are presented in [14] and in monograph form 

in mang. 

It should be noted that the results described above pertain to the interior and fully supported 
dissipation. The situation is much more delicate when the dissipation is geometrically constrained, 
where the essential support of the damping is localized to a subset of the spatial domain Q. In that 
case, the issue of propagating the damping from one area to another becomes the critical one. While 
this sort of problems has been previously studied in the context of stabilization to equilibria [MI f2~4"l |2"5] . 
the estimates needed for attractors are much more demanding. Previous methods developed in the 
context of stabilization no longer apply. Some long-time behavior results with boundary damping 
are presented in [131 114] , wherein nonlinear dissipation on the boundary acting via free boundary 
conditions is considered. These works, however, impose the rather stringent geometric restrictions of 
the entire boundary being star-shaped. Such restrictions are removed in [12] . where dissipation via 
hinged boundary conditions is considered; however this is done at the expense of limiting the class of 
dissipation to those of linearly bounded type. This restriction is needed since the elimination of the 
geometric condition is achieved via microlocal estimates [32] . which in turn force velocity dependent 
nonlinear terms to be linearly bounded. 

This brings us to the main contribution of the present manuscript. Our goal is to show that 
the fully nonlinear damping with essential support in an arbitrarily small layer near the boundary 
provides not only the existence of compact attractors but also desirable properties such as smoothness 
and finite dimensionality. Thus the original hyperbolic-like non-smooth flow is asymptotically reduced 
to smooth and finite dimensional dynamics. The result is valid for all types of boundary conditions 
with geometrically constrained dissipation, which can be nonlinear of any polynomial growth at infinity 
and with no restriction on the size of the damping parameter. 

We obtain this result by proving that the dynamics are quasi-stable - a concept introduced in [10] 
and [T2]. The ability to show quasi-stability is dependent upon: (a) a new method of localization of 
multipliers that allows smooth propagation of the damping from the boundary collar into the interior 
(even in the presence of boundary conditions (free) that do not comply with the Lopatinski conditions 
|4"2"] ) and (b) "backward" smoothness of trajectories from the attractor - a method used also in [2] 
and in [151 [5] - the latter in the context of geometrically constrained dissipation for wave dynamics. 

Lastly, we would like to note that while some of the methods developed for boundary dissipation 
[T2l [T4"] can also be used in the case of partially localized dissipation and Dirichlet - clamped boundary 
conditions, this is not the case with Neumann type (free) boundary conditions which violate strong 
Lopatinski [42] condition. In this latter case, propagation of the damping from the boundary layer via 
boundary damping estimates is obstructed by the well known lack of sufficient regularity (the absence 
of so called "hidden" regularity [3T] ) of boundary traces corresponding to the linear model [32] . Our 
aim in this paper is to develop a method which is effective for all kind of boundary conditions and does 
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not depend on hidden regularity, where the latter restricts the analysis to Lopatinski type of models. 
The key element for this are suitably localized multipliers estimates. 



1.3 Statement of Results 



Equipped with well-posedness of finite energy and regular solutions corresponding to ( 1.1 ) under one of 
the boundary conditions (C), (H), or (F), we are now ready to state our main results pertaining to long 
time behavior of solutions. In order to do this, we shall introduce the following unique continuation 
condition, denoted U C. We say that the system satisfies the U C property iff the following implication 



is valid for any weak solution (u,Ut) to (1.1): There exists T > such that 



u t = a.e in supp d x (0, T) =>■ u t = a.e in f2 x (0, T). 
It is clear that the UC property holds if d(x) > a.e. in fl. 

Remark 1.3. It is worth noting that due to the non-local character of von Karman nonlinearity, the 
unique continuation property for the von Karman plate is not fully understood. A now classical set 
of tools developed for plate equations and based on Carleman estimates [HI El [23 EH] do not apply. 
The non-locality of the von Karman bracket prevents propagation across the entire domain of weak 
damping localized to a small set. Therefore, we have the question: if the damping in the equation 
(represented by d(x)u t ) is zero in an open set of positive measure inside ofQ, does this imply that the 
solution u must also be inQ?; it remains open. In relation to our analysis here, if the general unique 
continuation property holds for the von Karman plate, then it immediately strengthens our result by 
allowing d(x) to vanish away from an open collar of the boundary. However, at present, the best we 
can state is a sufficient condition, namely that c?(x) > a.e. in order to satisfy the UC property. 

In addition, we shall assume validity of an asymptotic growth condition from below imposed on g(s). 
Such condition is typical |34j and necessary in order to obtain uniform decay rates of solutions in 
hyperbolic-like dynamics. It allows control of the kinetic energy for large frequencies. 

Assumption 1. There exist positive constants < m < M < oo and a constant p > 1 such that 

m < g'(s) < M\s\ p , \s\ > 1 
We now state the primary result in this treatment. 

Theorem 1.3. Take Assumption [7] to be in force. Let supp dDu and d(x) > «o > in uj, where 
lo CC O is any full collar near the boundary T. Then for all generalized solutions corresponding to 
solutions with initial data ||(mo, < R> there exist compact attractor A# £ H. If in addition, the 

UC property holds, then said attractor is global, i.e A^ = A for all R > 0. 

Remark 1.4. The UC property is needed in order to construct a strict Lyapunov function for the plate. 
In the absence of this Lyapunov function, the obtained results are of local character - as in the first 



part of Theorem 1.3 The UC property allows us to conclude that local attractors coincide with a 
global one. 

Theorem 1.4. In addition to Assumption 1 and the UC property, assume that there exists m, M > 0, 
and 7 < 1 such that < m < g'(s) < M[l + sg(s)]' y , for all seM. Then, 
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(a) the attractor A is regular, which is to say A C H (fl) x H 2 (fl) is a bounded set in that topology. 

(b) The fractal dimension of A. is finite. 

IP 

Remark 1.5. If we consider g(s) — |s| p+1 , then one can show that 7 = satisfies the above 

p -\- 2 

condition. 

There arc three main difficulties/novelties pertaining to the proof of the results stated above: 

(a) The nonlinear source is of critical exponent (lack of compactness). 

(b) The essential damping is geometrically constrained to a small subset lo. 

(c) The damping is genuinely non-linear (any polynomial growth at the infinity is allowed). 
These three difficulties are well-recognized in the context of studying long time behavior of hyperbolic- 
like systems where there is no inherent smoothing mechanism present in the model. In order to provide 
some perspective, it helps to add that geometrically constrained damping forces to use higher order 
multipliers which become supercritical when dealing with energy terms and nonlinear critical terms. 
Thus, any successful approach must rely on suitable cancellations, which must be uncovered for the 
specific dynamics in question. 

Similar issues appear when dealing with nonlinear damping. The damping term must be critical 
(in hyperbolic dynamics) in order to be effective (we recall that the essential spectrum of an operator 
can not be altered by a compact perturbation). The property of monotonicity of the problem does 
help when dealing with a single solution at the energy level. However, when dealing with long-time 
behavior, the protagonist is not a single solution but the difference of two solutions. In the study 
of the corresponding dynamics at the non-energetic levels (resulting from multipliers), monotonicity 
is destroyed. There is a "spillover" of the noncompact (in fact, supercritical) damping that must be 
absorbed. For this issue, different mechanisms need to be discovered (e.g. backward smoothness of 
trajectories, compensated compactness, etc). 

While recent developments in the field provide tools enabling us to handle a combination of any two 
of the difficulties listed above, the inclusion of the third prevents us from utilizing existing mathematical 
technology. The principal contribution of this treatment is to develop method which is capable of 
dealing with all three aforementioned difficulties simultaneously. The main ingredients of this new 
approach are (i) a localization method which allows us to show propagation of the damping without 
any requiring that the Lopatinski condition be satisfied, and (ii) backward smoothness of trajectories 
from the attractor with geometrically localized dissipation. 

We conclude this section by listing few problems that are of interest to pursue and still open. 

(1) C°° smoothness of attractors in the presence of nonlinear damping 

Regarding the first item, C°° smoothness of an attractor can be proved under certain restrictions on 
nonlinear damping by methods developed in |12j and also in an influential paper |23j . The treatment 
of a fully nonlinear and monotone dissipation is still not fully understood. 

(2) Damping restricted to a portion of an open collar 

Secondly, dissipation localized to part of the collar could be considered by assuming certain geo- 
metric conditions imposed on the uncontrolled part of the collar. Certain ideas presented in [T51 [5] 
should prove useful. 

(3) The UC property for a larger class of dampings 

Lastly, the third item is an wide open and important problem. Other mild forms of dissipation 
- such as viscoelastic weak damping - are also valid alternatives, however, a full understanding of 
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this problem would require the development of a substitute for Carleman estimates, which are not 
applicable due to localized structure of the von Karman bracket. The associated problem is related to 
controlling low frequencies - an endemic problem principally associated with strong stability [301 129] . 
It is worth noting that this problem is non-existent in the case of other nonlinear plates, such as 
semilinear plates with local nonlinearities or even Berger's model [6], where some forms of Carleman's 
estimates do apply [TJ [TH] . 



2 Long-time Behavior of Dynamical Systems 

In this manuscript we will make ample use dynamical systems terminology (see [21 [38l |4TJ [3 [T2] ) ; let 
(H,St) be a dynamical system with Jf = {x £ % : S t x = x for all t > 0} the set of its stationary 
points. 

We say that a dynamical system is asymptotically compact if there exists a compact set K which 
is uniformly attracting: for any bounded set D c H we have that 

lira d H {S t D\K} = (2.1) 

t— >+oo 

in the sense of the Hausdorff semidistance. 

(H,S t ) is said to be asymptotically smooth if for any bounded, forward invariant (t > 0) set D 



there exists a compact set K C D such that (2.1 1 holds. An asymptotically smooth dynamical system 
should be thought of as one which possesses local attractors, i.e. for a given ball Br(x) of radius R in 
the space T-L there exists a compact attracting set in the closure of Br(x), however, this set need not 
be uniform with respect to R or x £ H. 

A global attractor A is a closed, bounded set in % which is invariant (i.e. St A = A for all t > 0) 
and uniformly attracting (as defined above). 

A strict Lyapunov function for (H,St) is a functional $ on H such that (a) the map t — > <S>(S t x) 
is nonincreasing for all x £ H, and (b) $(Stx) — $(a;) for all t > and x £ W implies that x is a 
stationary point of (H, St). If the dynamical system has a strict Lyapunov function, then we say that 
(H,St) is gradient. 

In the context of this paper we will use a few keys theorems (which we now formally state) to prove 
the existence of the attractor. (For proofs and references, see [12] and references therein.) First, we 
address attractors for gradient systems and characterize the attracting set: 

Theorem 2.1. Suppose that (H,St) is a gradient, asymptotically smooth dynamical system. Suppose 
its Lyapunov function <fr(x) is bounded from above on any bounded subset of T-L and the set = {x £ 
% : <f>(x) < R} is bounded for every R. If the set of stationary points for (H,St) is bounded, then 
{%, S t ) possesses a compact global attractor A which coincides with the unstable manifold, i.e. 

A = Jf u {jV) = {x £ H : 3 U(t) £%, V t 6l such that 17(0) = x and lim d H (U(t)\J^) = 0}. 

t— > — oo 

Secondly, we state a useful criterion (inspired by [27] ) which reduces asymptotic smoothness to 
finding a suitable functional on the state space with a compensated compactness condition: 
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Theorem 2.2. Let (7i,S(t)) be a dynamical system, 7i Banach with norm \ \ ■ \ \. Assume that for any 
bounded positively invariant set B C 7i and for all e > there exists a T = T Ci b such that 

\\StXi - S T X 2 \\-H < ^ + ^e,B,T(xi,X 2 ), X { € B 

with a functional defined on B x B depending on e, T, and B such that 

lim inf lim inf ^ e r b (x m ,x n ) = 

m n 

for every sequence {x n } C B. Then (71, St) is an asymptotically smooth dynamical system. 

In order to establish both smoothness of the attractor and finite dimensionality a stronger estimate 
on the difference of two flows is needed. 

Theorem 2.3. Let x\, x 2 € B C 7i where B is a forward invariant set for the flow StXi. Assume that 
the following inequality holds for all t > with positive constants C\(B), C 2 (B),ujb 

\\S t x 1 - S t x 2 \\ 2 n < dWe-^Wx! - x 2 \\l + C 2 (B) max \\S tXi - S T x 2 \\ 2 Hi (2.2) 

rS[0,t] 

where 7i C 7i\ is compactly embedded. Then the attractor A associated with the flow St posesses the 
following properties: 

(a) The fractal dimension of A. is finite. 

(b) For any x £ A one has ~^{Stx) € C(R,7i). 

Remark 2.1. The estimate in ( |2.2[ ) is often referred to as a "quasistability" estimate. It reflects the 
fact that the flow can be stabilized exponentially to a compact set. Alternatively, we might say that 
the flow is exponentially stable, modulo a compact perturbation (lower order terms). We assert that 
the lower order terms being quadratic is important for the validity of Theorem |2.3| 



The proof of Theorem 2.3 employs the idea of "piecewise" trajectories introduced in [Sol 40 . This 
allows to generalize previous criteria for finite-dimensionality [51 1441 1411 117j by reducing the problem 
to validity of quasistable estimate. 

2.1 Approach and Outline of the Paper 



To show our main result on the existence of the global attractor for (1.1 1 with boundary conditions 



(C), (H), or (F) we make use of the theorems above. First, we note that in the case of any boundary 



conditions, the von Karman system in (1.1) is gradient with Lyapunov function S(t) (under the 
assumption that the UC property is satisfied, e.g. in the case that d(x) > a.e. in f2). We refer 
to [12] for the details. Moreover, the set of stationary points for the dynamical system generated by 



<fTTT])(C), ([rTJ>(H), or ([13]) (F) is bounded. This latter fact follows from (see Q2]). Hence we are 



in a position to use Theorem |2.1| if we can obtain an inequality of the form in Theorem |2.2| to show 
asymptotic smoothness of the system. We will analyze z, taken to be the difference of strong solutions, 
and make use of the linear energy E z (t) — ||Az|| 2 + ||z t || 2 (then, via a standard limiting procedure 
obtain our estimate for generalized solutions as well); this estimation will produce our functional ^ in 
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Theorem 2.2 Our main tool in estimating E z (t) will be the use of two multipliers: f\Z and /i-V((/2)z)), 
where h will be a suitably chosen C 2 vector field and ft are appropriate localization functions. 

First, we perform multiplier analysis as generally as possible, without imposing boundary condi- 
tions. Later on, we shall use boundary conditions (either clamped, or hinged or free) in order to obtain 
the smoothness inequality in Theorem |2.2| 

After establishing the existence of the attractor, we proceed to show that it has additional regularity 
than that of the state space, and also that it has finite fractal dimension. The ultimate goal is to 
prove a "quasistability" estimate for the difference of general trajectories and apply abstract Theorem 
|2.3[ however doing so directly in this case is difficult. Instead, we will prove a modified quasistability 
estimate (via similar methods in the asymptotic smoothness calculation) which applies to the difference 
of two terms along the same trajectory (i.e. the difference of u(t+h) — u(t)). This term is substantially 
easier to analyze, since we have continuity in t in T-L and both u(t + h) and u(t) converge to the same 
point of equilibrium for t — > oo and t —¥ — oo. Our modified quasistability estimate hinges upon the 
bounding terms being quadratic, so upon division by h and taking h \, we can show additional 
regularity of elements from the attractor for sufficiently negative times T << 0. Proving this will 
depend upon a trajectory being 'close' to a point of equilibrium, and hence yielding 'smallness' of 
the velocity of the solution. Propagating this regularity forward via the dynamical systems property 
will then allow us to show the additional regularity of the attractor. We will then proceed in a 
standard fashion to show that this additional regularity of the attractor yields the true and sought 



after quasistability estimate in Theorem 2.3 which will produce the finite fractal dimension of the 
attractor. 



3 Asymptotic Smoothness 



In this section we prove that the dynamical system generated by (1.1 1 is asymptotically smooth. We 
will refrain from imposing boundary conditions until absolutely necessary in the hope of unifying the 
treatment of (C), (H), and (F). 

Note that the new variable z = u — w, where (u(t),Ut(t)) = St(uo,ui), (w(t),wt{t)) — St(wo,W\) 



are solutions to (1.1 ) with initial data taken in bounded set in B C H. On the strength of Lemma 1.1 



and lip we may assume that there exists R > such that 

\\(u(t),u t (t)\\ n <R, \\(w(t),w t (t)\\ n <R, i>0 (3.1) 

The difference of two trajectories z = u — w solves the following PDE: 

z tt + A 2 z+g{z)+ F{z) =0 in Q, (3.2) 
z(0) = u -w ; z t (0) =U!-Wi 

where 

F(z) = -(fv(u) - fv(w)), and Q{z) = d(x)(g(u t ) - g(w t )) 

The above evolution is equipped with appropriate boundary conditions (C), (H), or (F) which will be 
specified later. 
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3.1 Multipliers 

Ultimately, we will need a pointwise bound (in time) on the functional E z (t) as defined above. To 
achieve this bound, we will employ multiplier methods based on specially chosen cut-off functions A 
and fi. These functions are taken to be C°°(f2). Later, we will choose the supports of the derivatives of 
A and y, to be contained in the damping region w, where the damping g(u t ) is effectively localized; the 
cut-off functions will be chosen in this way so as to reconstruct the full energy E z (t) via the multipliers, 
bounded in terms of the damping. However, for now, we can consider supp A C to be arbitrary. 

We define the variables <f> = Xz and ip = fxz. The use of the cut-off functions will produce 
commutators active in the regions of u> where the cut-off functions are non-constant. Lastly, we will 
make use of the following notational conventions. First, to describe (a) lower order terms: 

l.o.tJ = sup II/WIIL,, l-o-t.{ = sup ||/(t)||a-„, 

[0,T] ' [0,T] 

where < rj < 1/2, and (b) boundary terms: B.T. f = {Afd„f - <9„(A/)/j 

Remark 3.1. We note that the use of different notations for lower order terms is necessary in the 
handling of dissipation estimates. Specifically, we must treat the dissipation terms differently when 
dealing with asymptotic smoothness type estimates, and the estimates which will ultimately yield the 
quasistability estimate. 

3.1.1 Multiplier 

Let P and Q be two differential operators. We will make use of the commutator symbol given by 

[P,Q]f = P(Qf)-Q(Pf), 



We shall work with smooth solutions guaranteed by Theorem 1.2 Multiplying the PDE in (3.2) by A 
we arrive at 

4> tt + A 2 + XQ{z) + A-F(z) = [A 2 , \]z. 

Now, we employ the multiplier <j). This is an equipartition multiplier which allows us to reconstruct 
the difference between the potential and kinetic energies. The following Green's identities are available 
|34j for sufficiently smooth functions z and <j): 

/ A 2 zcf> — / AzA(f> + / (d v Az4> — Azd v (t>), clamped and hinged B.C 
Jq Jn Jr 

[ A 2 z<f> = a(z,(j>) + (3 [ z 3 <J)+ [ (B 2 Z(f)-Bizd v ^), free B.C 

.Jn Jr! Jr! 

Using the first formula for clamped or hinged boundary conditions yields: 
/ {|A0| 2 -|0 t | 2 } = f [A 2 ,A]z0- f \{g{z)+F{z)} ( j > + [ {A00„0-ft,(AM}-(&,0)|o (3.3) 
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Making use of standard splitting and Sobolev embeddings, we arrive at 

'' {||A () |! J -iNH J ) / B.T.-+ I l[A-.A]:)o + ^ 
In the case of free boundary conditions, the equipartition of energy takes the form 



/ T {||A0|| 2 -||0 t || 2 } < f B.T.*+ [ ([A 2 ,A]z)0 + / \{g{z)+F(z)}^ + C{E{T) + E{Q)) (3.4) 

JO J£ JQ JQ 



{a(0,0)+/?Hl 4(r) -|N| 2 } < / (Bi#-S 2 <R0) + / ([A 2 ,A]^)</» 

JSi JQ 

+ / \{g{z)+F{z)}<t> + C{E{T) + E{Q)) (3.5) 
JQ 

We note for all boundary conditions (C), (H), the boundary terms B.T.^ = 0. In the free case (F) we 
have B\4> — 0, IBify — 2f3<t>uw where the latter term contributes a lower order term to the estimate. 

To continue with our observability estimation, we must explicitly bound the commutator / [A 2 , \]z<j> 

JQ 

Purely algebraic calculations give 

[A 2 ,A]/=A 2 (A/)-AA 2 / 

= (A 2 A)/ + 2AAA/ + 2(VA, V(A/)) + 2(V(AA), V/) + 2A(VAV/) (3.6) 

The calculation above implies that the commutator [A 2 , A] is a differential operator of order three. In 
order to exploit this in the calculations with the energy, we need to reduce the order of differential 
operator acting on a solution via integration by parts. This is done below. 

This computation makes sole use of Green's theorem. For the sake of exposition, we do not impose 
any boundary conditions: 

VAu(0VA) = - / (Au)div(<£VA) + / (0Au)VA ■ v (3.7) 
a Jn Jr 

/ A(VAVw)0=- / V(VAVw)V0+ / d v {VuVA)(j> (3.8) 
Jn Jn Jr 

Note that due to the fact that the support of VA is away from the boundary, all of the boundary 



terms in the above expressions (|3.7|) and (|3.8|) will vanish. Moreover, 

VAVAw<?!> + / A(VAVu 



<Ca|H| 2 ||0||i (3.9) 



Hence to conclude our <f> multiplier estimate, we have the following technical lemma: 



Lemma 3.1 (Preliminary Estimate). Let <p = Xz, as defined above, where z solves ( 3.2) with boundary 
conditions (C) or (H). Then, there exists < C < oo such that 



r-T 

|2 1 1 i ||21 



/ {||A0|| 2 -||0 f || 2 } < C(T,X)l.o.t. z + \{Q{z)+F{z)}d? + C(E z (T)+EM) (3-10) 
lo Jq 
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In the free case (F) 



j {a{^4>) + p J^-W&W 2 } < C(T,\,R)l.o.t. z + J \{g(z)+j u (z)}(t> + C(E z (T)+E z (0)) (3.11) 



Proof. Taking into account (3.9 1 in (3.4), we have 



{||A^|| 2 -||0 t || 2 } < C(T,X)l-o.t. z + / \{G(z)+T(z)}0+ / BT^ 

n JQ JS 

+ J {d v {A(\z)){\z)-A(\z)d v {\z)-d v (Az)\ 2 z 
{Az)d v z + A 2 (Az)d„z} + C(E(T) + E(0)) 



2Xz( 



(3.12) 



Taking into consideration boundary conditions (C) or (H) in (3.4), noting that B.T^ — and ac 



counting for the fact that the boundary terms resulting from the commutators vanish leads to the 



first statement in the Lemma. Calculations in the free case are analogous, and result from (3.5) and 
Bi4> — 0, 82(f) — 2fi(j)uw, where the latter term contributes a lower order term to the estimate: 

£ 2 < 



< 2/3 f H 2 MM < 2/3R 2 \\(t>\\j < C(R)l.o.t z 



□ 



3.1.2 Multiplier 2: h ■ 

For the first part of this section, we specify only that supp fi D T — 0; otherwise, we keep fi as general 
as possible, specifying it at the last possible moment. Additionally, define a set M = supp V/i = 



{x 6 Q\ii ^ constant}. Now, if we multiply (3.2) by /i, and recall that i/j = /iz, we obtain 

tPu + A 2 V + plQ{z) + iiT{z) = [A 2 , fj]z 

where Q{z) = d(x) (g(u t ) — g(w t )) and ^{z) = —(fv(u) ~ fy(w)), as before. We now make use of the 
multiplier h ■ V?/>, which we write as h\7ip henceforth; there are various choices for the vector field h, 
situationally dependent, however here we need only take h = x - x £ M 2 in order to obtain control 
on the potential energy of the plate. Now, as in the previous section, we multiply the last equality by 
our multiplier and use Green's Theorem to obtain 

(|V t | 2 + |AVf) < C(E Z (T) + E z (0)) + f ^g{z)(hVi>)+ [ /i.F(z)(WV0+ / [A 2 , ^z(W^) 

JQ JQ JQ 

By explicitly writing out the commuator, and taking into account the support of V/i, upon splitting 
we obtain: 

[A 2 ,^]z(W^)= I I [A 2 ,^]z(WV) < C(aO I f |Az| 2 +C(T, \i)l.o.t. z (3.13) 

Q JQ JM JO JM 



Now, at this point we specify the specific structure of the supports for A and jx (which up to now have 
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been general). The following picture illustrates our choice for these supports and their relationship to 
the damping region u: 




We emphasize that (a) the set M C supp A and (b) supp A and supp n overlap inside the damping 
region ui and that supp A U supp(,u) = il. Since M c {x e ft : X(x) = 1}, we have the following 
inequality: 



K 2 



[A 2 ,ii]z(hVip) <C(fi) / \Az\ 2 + l.o.t 

Q JO JM 

<C([x) '' 



f I \Az\ 2 + C(n,T)l.o.t. z 
Jo J \=i 

<C(jjl) [ \A</)\ 2 + C(n,T)l.o.t. z (3.14) 



3.2 Energy Recovery Estimate 

We may now appeal to our calculation with the <f> multiplier previously, to obtain our preliminary ip 
estimate: 



Lemma 3.2 (Preliminary ijj Estimate). Let ip = fiz, as defined above, where z solves (3.2) with any 



boundary conditions under considerations. Moreover, assume supp (/i) is bounded away from V. Then, 
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in the case of clamped (C) or hinged (H) boundary conditions we have 

(|^| 2 + |A^| 2 ) < C{n,\){{E,{T) + EM) + f Q ^{ g ^ +F(*)}(hVtl>) 

+ [ \{Q(z)+F{z)}^ + C{T)l.o.t. z (3.15) 
In the case of free boundary conditions (F) 

+ a(V>,^))< C(n,\){(E z (T) + E z (0))+ f »{G(z) + F(z)}(hViP) 



1 



+ j \{Q(z)+F(z)}^ + C{T)l.o.t. z } (3.16) 



Q 

We note that in the nonlinear boundary term associated with the operator B2 vanishes due to the 
fact that the support of [i is away from the boundary. 

We may now combine the estimates from Lemma [3~T] and Lemma [372] to obtain an estimate on the 
total energy (with either form of boundary conditions (C) or (H) or (F)): 

{Ut\\ 2 + \\M 2 +^^)+<^,^)+P / H 4 }< C(n,X){(E z (T)+E z (0))+ f i x{Q{z)+F{z)}{hV^) 

JT t L JQ 

\{G(z)+J r (z)}4> + J \\^t\\ 2 + C(T,R)l.o.t. z } (3.17) 

By our choice of supports for fi and A we note that the LHS of the above equation overestimates 
the total energy E z (t). On the RHS of the estimate we have the term Jq \4>t\ 2 , which we replace by 

Jo JL l z *l 2 smce SU PP A C oj and on supp A, A < 1, so we have that 



\m 2 < r 1 w 2 



Q 



Making the appropriate changes above in Lemma 3.1 and Lemma |3.2| we have the analogous result 
for the free boundary conditions (F). Hence we can conclude 

Lemma 3.3 (Preliminary Energy Estimate). For any boundary condition (C), (H), or (F) we have 
E z (t)< C(n,\){(E z (T) + E z (0))+ f ti{g(z)+F(z)}(hVi>)+ f \{Q{z) + F{z)}(t> 



T 

2 



\z t \ l +C(T)l.o.t. z (3.18) 

Remark 3.2. At this point, we impose clamped (C) or hinged (H) boundary conditions, in order to 
simplify and streamline the analysis. At the end of this section, we discuss the boundary conditions 
(F). 
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If we take into account the supports of A and /j, (dropping dependence of the constants on fi, A, and 
Q) then (3.18) with clamped boundary conditions becomes 



/ E z (t)< c\e z (T)+E z (0)+ f {g(z)+F(z)}{hW^)+ [ {g{z)+F{z)}z\ 
Jo L Jo Jo ' 



C \ \ \z t \ 2 + C(T)l.o.t. z (3.19) 



Remark 3.3. At this point we pause to point out that the estimate we have shown above in (3.191 
will be used in the sections to follow, specifically in the quasistability estimate. In particular, we must 
handle the damping terms (involving ut, Wt) differently in the estimation for asymptotic smoothness, 
versus the estimation for quasistability. 

By the assumptions on g in Assumption [l] for every S there exists Cs > such that 

\u t - w t \ 2 < 5 + C s (g(u t ) - g{w t )) (u t - w t ) . 

This gives that 



\zt\ 2 < TS\n\ + C{5) / / (g(u t ) - g(w t ))z t 

Jo Jui 

or, simplifying, and taking into account lj C supp d and that d(x) > oq > 0, we have 

\z t \ 2 <5 + c(s,T,n) f g(z)z t 







So taking into account the last inequality in (3.19) , we obtain 

l-T 



E z (t)<6 + C{E Z (T)+E z (0) + j o j (G(z)+T(z))z 

+ f (g(z)+?(z))hVii> + c(5,T) ( g{z)z t X 

Jq Jq j 

+ C(T)l.o.t. z 



where the constant C does not depend on T. Recall, u and w are solutions to (1.1) corresponding 
to some initial conditions y\ and y2, satisfying Styi = (u(t),ut(t)) and Sty% = (w{t),wt{t)) for the 
evolution St associated to the plate dynamics. We can assume that yi, 1/2 6 Wr for some R > R*, 
where the invariant set Wr = {(u, v) g H, §{u, v) < R} . Assuming the solutions u and w are strong, 
by the invariance of Wr we have 

||u(t)|| 2 + |Mt)ll + lk(*)ll2 + IM*)||<C(fl), t>0 (3.20) 

\\u(t)\\c(n) + \Ht)\\ c{ n)<C(R), t>0 (3.21) 

Recent developments in the area of Hardy- Lizorkin spaces and compensated compactness methods 
allow one to show the following 'sharp' regularity of the Airy stress function v. 
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Theorem 3.4 (Sharp Regularity of the Airy Stress Function). 

\\v(u)\\ W 2. x < C\\u\\l, \\v(u,w)\\ W 2,o* < C||w|| 2 |H| 2 

where we have denoted v(u,w) = — A~ 2 [u, w] and £>(A 2 ) = H 4 (fl) n Hq({Y) 
Making use of the above inequalities, we have the estimate 

\\[v(u),z}\\<C\\u(t)\\l\\z\\ 2 <C(R)\\z\\. 

Additionally, we have 

\\v(u) - v(w)\\ W 2,oa = \\v(z,u + w)\\ W 2,co < C\\z\\ 2 (\\u\\ 2 + \ \w\\ 2 ) ■ 

Therefore, 

||J-(z)|| = \\[u,v(u)}-[w,v(w)] + [z,F ]\\ = \\[v(u)-v(w),z] + [v(w),z] + [z,F ]\\ < C(R)\\z\\ 2 , t > 0. 
So we obtain 

' F(z)z < ( T(z)z < e { \\z(t)\\ 2 2 dt + C{T,e)l.o.t. z (3.22) 



and similarly 

r T{z)hVij<C{R) I ||*(t)|| 2 ||V(i)lli < * / Mt)\\ 2 2 + C(T,e)l.o.t. z , (3.23) 



Q 

(where again, dependence of constants on Q,,lj, and h are supressed). To proceed, we need estimates 
on the dissipation. By the energy equality 

E Z (T)+ [ [ G(z)z t =E z ( S )+ [ [ F{z)z u (3.24) 
Js Jo. Js Jn 

we have 



Q(z)zt<C(R)+ / T{z)z t (3.25) 

JQ 

Taking into account the embedding i? 2 ~''(£l) C C(f2) for < r\ < 1, we see 

G(z)z< [ d( X )(\g(u t )\ + \g(w t )\)\z\ 
> JQ 

< C\\z\\ c(0 . T , c(n)) [ d(x)(\g(u t )\ + \g(w t )\) 



< C\\z\\c(o,T;H 2 -i(n)) / d(x) (\g(u t )\ + \g(w t )\) . 

JQ 

Splitting the region of integration according to \ut\ < 1 and \ut \ > 1, and similarly according to \w t \ < 1 
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and \w t \ > 1, we obtain 

d(x)(\g(u t )\ + (g(w t )\) < g(l)\\d\\ Loo{n) meas(Q) + / d(x) (g{u t )u t + g(w t )w t ) <C{R,T) 



Hence 

/ g(z)z<C(R,T)l.o.t. z 1 (3.26) 
JQ 

Now applying Holder's inequality with the exponent r > 1 we see 

e(*)wv><csu P ||vv(i)||r' / ^( x ) r (l5K)l r + l5K)l r ) 

[0,T] JQ 

where — I — ^ = 1. Taking r = 1 H , and again splitting the region of integration according to 

r r p + 1 

|itj| < 1 and |ut| > 1, and using the polynomial growth condition imposed on g in Assumption 1, we 
obtain 

' d( X ) r \g(u t )\ r < C(d){ 5 (l)meas(Q)+ / d{*)g{u t )u t } < C(R)(T + 1) 

' JQ 

Since the same computations hold for terms in w, and we have the continuous embedding H 1 ^ (fi) 
L r > (ft) for sufficiently small S, we have 

/ g(z)KVil) <C{R i T)l.o.t.{ (3.27) 
JQ 

Hence by the above estimates, we have 

E z (t)<C{E z (T)+E z {0) + S + C{R,5) + C(S) f F{z)z t + C{R,T)(l.o.t. z + l.o.t.f)} 
o Jq 

and eventually by (3.24) we have 



T 



E. 



ii 



'»(t) < C t i [ E{T) + 8 + C{R,8)+C{8)\ j J T{z)z t | + C(R, T)(l.o.t* + l.o.t.\) } (3.28) 



where we write C* to emphasize that this constant does not depend on T. If we integrate (3.24) over 
(0, T) with respect to the variable s, and take into account (3.28), we may choose T sufficiently large 
(T > 2C*) and e sufficiently small (with respect to T) such that 

Lemma 3.5 (Asymptotic Smoothness Estimate). 
E Z {T) < 6 + M + | jf J-(z)z t | + \£ £ J n F(z)zt\J + C(e, R, T)(l.o.t. z + l.o.t.{) (3.29) 

We are now in a position to prove Theorem |1.3| on the existence of a compact attracting set A. 
For this we shall invoke the abstract Theorem 12.21 

Completion of the proof of Theorem [7 
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To apply Theorem 2.2 we need to construct a functional $ £ ,.r,t such that 

liminfliminf$ eifljT (2/ rl ,y m ) = 



m— too n— too 



for every sequence {y n } from B (following from Theorem 2.2 1. The functional will contain "noncompact 
and not-small" terms in the inequality (3.29). More specifically, for any initial data Uq — (uq, iti), Wq = 
(wq, Wx) € B we define 



$e <R ,T(Uo,Wo) = / (f{z),Zt) 



T r T 



where the trajectory z = u — w has initial data Uq — Wq. The key to compensated compactness is the 
following representation for the bracket: 

(F(z),z t ) = ~{ - \\Av(u)\\ 2 - \\Av(w)\\ 2 + 2([z,z],F Q )} - ([v(w),w],u t ) - ([v(u),u],w t ) (3.30) 

Integrating the above expression in time and evaluating on the difference of two solutions z n ' m = 
w n — w m , where w' — w yields: 



n— too m— >oo 



lim lim / (J-(^-),^" l )=-{||A W H(t)|| 2 -||A«H(T)|| 2 } 



(3.31) 



lim lim / {([v(w n ) 1 w n },w] n ) + ([v(w m ),w m ],w?)\, 



where we have used (a) the weak convergence in H 2 (£l) of z rhm to 0, and (b) compactness of Av(w) 
from H 2 (Q) — > Li2(£i). The iterated limit in (3.31) is handled via iterated weak convergence, as follows: 

lim lim [ T {([v(w n ),w n ],w?) + ([v(w m ),w m }, W ?)} 

= 2j\[v{w),w]^t) = \\\Av{w)(t)\\ 2 - l -\\Av{w){T)\\ 2 . 
This yields the desired conclusion, that 

lim lim / (F(z n > m ),z?' m ) = 0. 

n— ¥oo m— ¥oo J ^ 

The second integral term in $ is handled similarly. As a consequence we obtain 

liminfliminf$ ei fl iT (yn,J/m) = 0. 

m— >oo n— >oo 

Now, we define 

$e,R,T = $ + {l.0.t. Z +l.0.t.l), 

and noting that the terms (l.o.t. z + l.o.t.f ) in (3.29) are compact with respect to H 2 (£l) via the Sobolev 



21 



embeddings, the final conclusion follows by taking T sufficiently large and e sufficiently small. This 
concludes the proof of smoothness estimate required by Theorem |2.2| Thus, the dynamical system is 
asymptotically smooth. In addition, stationary solutions are bounded (due to the inequality in Lemma 
1.1) and the set {(uo,Ux) £ H, <?(uo,ui) < R} is positively invariant. This implies ([H]) the existence 
of local attractors Ar - the first statement in Theorem |1.3[ For the second statement we just note 
that the UC property renders the dynamical system gradient with a strict Lyapunov function. Thus, 
asymptotic smoothness and boundedness of the set of stationary solutions implies the existence of a 
compact global attractor. 



4 Regularity and finite dimensionality of the Attractor 

Let A be the global attractor corresponding to the flow St, as established in Section 3. To prove 
finiteness of the fractal dimension of the set A, we shall use Theorem |2 . 3| which is based on a quasista- 
bility estimate. It will turn out that the quasistability estimate will follow from a more direct estimate 
written for the difference of two solutions z = u — w at sufficiently negative times. Since the system 
is gradient, the latter due to the UC property assumed, we have that trajectories from the attractor 
stabilize asmyptotically to equilibria points for both positive and negative times. This is to say: any 
x £ A belongs to some full trajectory 7 = {(u(t),Ut[t)),t € R} and for any 7 C A there exists a pair 
{e, e*}cJV (set of equilibria) such that 

(u(t),u t (t)) (e,0) as t ->• 00, (u(t),u t (t)) -t (e*,0) as t -t -00 (4.1) 



4.1 Quasistability Estimate 

We shall follow a general program developed in [T3] and supported by PDE estimates derived in 
previous sections and specific to localized dissipation. 

With the previous notation, we state the following lemma which gives a preliminary estimate for 
quasistability inequality: 

Lemma 4.1. Let z = u — w where (u(t),ut(t)), (w(t),wt(t)) £ A with w(t) — u(t + h), < h < 1. 
Then, there exists T u £ K such that for all —00 < s < s + Tq < T u the following inequality holds: 



E z (s + T )+ T ° E z (t) < C(A)T Di +T » +C(A)T sup \\z(r)\\l_ v (4.2) 

Js tG[s,s+T ] 

for rj > 0, where 



D\l = / / d( X )(g(u t ) - g(w t ))z t 
Jt-i Jn 

This lemma will provide a tool for proving additional smoothness of the attractor. In fact, we will 



select a time T which will be negative (in line with the convergence in (4.1 \). The above estimate, a 
posteriori, will yield the quasistability estimate in Theorem |2.3[ and hence finite dimensionality of the 
attractor. More specifically, we shall follow the program described below after completing the proof of 



the above lemma: (1) Apply Lemma 4.1 with some negative time s + T, to conclude (by homogenity) 



that the solution at T has additional derivatives - and is in fact a strong solution. For this step, we will 
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exploit closedness of velocities to a point of equilibria, along with a decomposition of the von Karman 
bracket, which allows us to take the advantage of this closedness. (2) The above gain of regularity is 
then propagated forward on the strength of regularity theory for von Karman evolutions. This gives 
an algebraic embedding for the attractor, but boundedness of the attractor in the higher norm will 
remain to be established; (3) This boundedness, that is, the topological embedding, will follow from a 
covering theorem, and the compactness of the attractor (shown in the earlier section). These details 
are given in |12j . Thus, at the end of this procedure we can extend validity of the inequality in Lemma 
4.1| to all times — oo < s < t < oo, and we obtain the conclusion of the additional regularity of the 



attractor. This then allows us to infer the quasistability estimate in Theorem 2.3 (which is the same 



estimate as in Lemma 4.1 but valid for all times s < t and all pairs of trajectories (u(t),ut(t)) and 



(w(i),wt(t)) lying in the attractor. 

Thus the crux of the proof of regularity and dimensionality of the attractor reduces to the demon- 



stration of Lemma 4.1 We also note that in comparison with the asymptotic smoothness inequality, 
the inequality in Lemma |4.1| is more demanding. This is due to necessity of keeping at least quadratic 
forms in the lower order terms. This very demand forces the damping to have at least linear growth 
at the origin g'(0) > 0. (Such restriction is typical - if not necessary - whenever regularity or finite 
dimensionality of attractors becomes a concern.) 



4.2 Preparation for the Proof of Lemma 4.1 



In the proof of Lemma 4.1 we will make use of the recovery estimate in (3.19) and the energy relation 



(3.24) as our main tools. Beginning with (3.19), and taking into account estimates involving F{z) in 
(3.22), (3.23) and the linear growth condition g'(0) > in (3.19), we arrive at 



o Jn 



\z t ? < CDl(z) 



Applying the above inequality in (3.19) gives: 

i-T 



E- 



(t) < c[dI{z) + E Z {T) + £7,(0) + 



G(z)z 



+ 



G(z)hVz 



Q 



X + C(R,T)l.o.: 



(4.3) 



Now, in tackling quadratic dependence of the dissipation terms, we give the following proposition 

Proposition 4.2. Let assumptions of Theorem 1.4 be satisfied, and take z be a solution to (3.2). Then 
there exists S > such that 



Q{z)z 



< 5D%(z) + C(<5, i?,T) sup < rj < 2 - 7 

[0,T] 



(4.4) 



G(z)hVi 



<5D T (z)+C(5,R,T)sup\\z\\ 2 2 „ v , < v < 1 - 7 

[0,T] 



(4.5) 
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where G(z) = d(x) (g(u t ) — g{wt)) and Dq(z) = / Q(z)z t . 

Jq 

Proof. We note that the assumptions on the damping function g (namely, montonicity and the poly- 
nomial growth condition in Assumption 1) imply that 



g(s 2 ) -g(si) 
32 - si 

Using the Jensen inequality we estimate 



< C[l + g( Sl ) Sl +g(s2)s 2 r. 



(4.6) 



\z\<S\z t \+C(5) 



The above, along with (4.6), gives 



Q{z)z 



< 5D< (z) + C(S, M) / d(x) (1 + (g (u t )u t y + (g (w t )w t y) \z\ 



~f\ ui 2 



Now, applying the Holder inequality with exponent p — 7 1 and Sobolev's embedding H 2 r ' '(O) C 



L 2 (O), and taking into account energy equality (1.13) we arrive at 

1-7 



Q{z)z 



Q 



< 5D^(z) + C{S,R)l.o.t. z 



The inequality in (4.5) can be shown analogously. 



□ 



So, taking into account (4.4) and (4.5) in (4.3) we obtain 



E 



(t) < c[Dl{z) + E Z (T) + E z (0)} + C(R,T)l.o.t. 



We note that C does not depent on T, and l.o.t." is of quadratic order. By using semigroup property 
and reiterating the same argument on the intervals [s, s + T] one obtains 



T+s 



E. 



(t) < c{dJ +s {z) + E Z (T + s) + E z (s)} + C(R,T)l.o.t. z (s,T + s) (4.7) 



where l.o.t. z (s, T + s) denote lower order terms collected over the interval [s, T + s]. 

In order to prove (4.2), we have to handle the non-compact term (^(z), z t ). A technical calculation 
based on the symmetry properties of von Karman bracket gives us the following proposition whose 
proof is given in |12j . 



Proposition 4.3. If it, w G C([0, t]\ H 2 (il)) n C^QO, t]\ L 2 (ty) andz = u-w then 
where 



(4.8) 



Q(z) = (v(u) + v(w), [z, z\) — ||A«(u + w, z) 
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P(z) = - (u u [u, v(z, z)]) - (w t , [w, v(z, z)]) - (u t + w t , [z, v(u + w, z)]) . 
Now, we can state the following lemma: 
Lemma 4.4. Let u(t) and w(t) be two functions from the class 

C([M];ffo(^))nC 1 ([ s ,t];L 2 (0)) 

for s, t € K, s < t, such that 

\Hr)\\l + \\u t (r)\\ 2 < R\ \\w{r)\\l + |K(r)|| 2 < R\ re [s,t] 
Let z(t) = u{t) — w(t) . Then for n > 



[ Wz),*) <C(R) sup \\z(r)\\ 2 2 _ n +C(R) f 

Js T£[s,t] Js 



\m + \m\\\\z\\2- 



(4.9) 



(4.10) 



Proof. The inequality follows from the basic properties of von Karman bracket |12j and the decompo- 
sition in Proposition 4.3. □ 



4.3 Completion of the Proof of Lemma 4.1 



Let 7 = {(u(t), u t (t)) : t G M.} be a trajectory from the attractor A, and let < ft- < 1. It is clear that 
for the pair w(t) = u(t + h) and u(t) satisfy the hypotheses of Lemma 4.4 for every interval [s, t]. We 
shall estimate the energy E z (t) of z(t) = z (t) — u(t + h) — u(t). Here we critically use the estimates 
for the noncompact term involving J~{z). By (4.10), we have 

(J-(A^) <C(R) sup \\z h (T)\\j_ r) + C(R) [\\\u t (r + h)\\ + \\u t (T)\\)\\z h (r)\\l (4.11) 

r£[s,t] Js 

for all — oo < s < t < +oo. Since we have the characterization A = AI u (Af), where Af is the set of 
equilibria, we have 

tondHSWxivwiStWliV ) - for any W E H*(n) x L 2 (Sl) 
Hence, for any e > 0, there exists (independent of h but depending on the trajectory 7) such that 

IK(r)|| + ||ut(r + /i)|| ^elCiR)}- 1 for any t < T^. 
Taking into account the last inequality in (4.11), we arrive at 



(4.12) 



(F(z h ),z?) <C(R) sup \\z h (T)\\l_ v + e \\z h {r)\\l 

TE[s,t] Js 

for all —00 < s <t < T^. Using the energy relation (3.24), we find from (4.12) that 

E z (s)< E z (t)+ f f g(z)z t + C(R) sup ||2 h (T)||!_„ + e / ||« h (r)||ldr. 

Js JO, re[s,t] Js 



(4.13) 
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and similarly 



E,(t) < E z (s) + C(R) sup \\z h (T) 

re[a,t] 



I2— jj 



(4.14) 



for all —00 < s < t < T*. Now, if we apply (4.7) on each subinterval [s, s + To], we have 

"S+T 



(t)<c{d s + t °(z) + (E z ( S + T )+E z ( S ))} + C(R,T ) sup ||z(t)||1_ 

L ' re[s.s+T ] 



Taking into account (4.13) in the last inequality and choosing e sufficiently small we arrive at 

ps+T 

/ E z (t) <c{d s s +t <>(z)+E z ( S + T q ) + C(R,T ) sup ||z(r)||tJ 

Js L t£[s,s+T ] ' 

Now, integrating (4.14) and considering the previous inequality we have 

^ E z (t) < C{A,T )D S + T «{z)+C{A 1 T ) 

re[s,s+T ] 



E z (s + T Q )+ r +T ° E z (t) <C(A,T )Dl +T °(z)+C(A,T ) sup ||«(t)||1_, 

Js tG[s,s+T ] 



which gives (4.2) and thus proves Lemma 4.1 



5 Proof of Theorem 1.4 



5.1 Proof of Part (i) in Theorem 



1.4 



Having established Lemma |4.1| we now proceed with the proof of improved regularity for the attractor. 
This is done, as in [14] ■ in three steps: 

Step 1: Smoothness for negative times 

By (4.2) and energy relation (3.24) written on the interval [s, s+T ] we can choose a constant < fx < 1 
such that u h (t) = hr 1 z h (t) satisfies the following estimate 



E uk (s + T ) < nE uh (s) + C(T ) sup llu^a + T)!!!., 

re[0,T o ] 



(5.1) 



for all s < T 7 — To, where T 7 = (depending on the trajectory, but not h) for some eo > and 
Tq > 0. Now using interpolation, and taking the supremum over the interval (-co, T 7 — T ) we obtain 



sup E uh (t) < 

r£(-oo,T 7 ] 



1+fl 



sup E u h(t)+C(T ) 

r£(-oo,T T ] 



This implies that 

E u k (s) < C(T ) for all s e (-oo, T 7 ] 
After passing to the limit h->0we get 



(5.2) 



lh t (t)|| 2 + |Mt)||^<C for all se (-oo,T 7 



(5.3) 
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By (5.3), for u t £ H 2 (Q) C C(fi), we have g(u t ) £ (7(0) c L 2 (0) by the continuity of g. Hence, elliptic 
regularity theory for A 2 u — —u t t — d(x)g(u t ) with the boundary conditions gives that ||u(i)||4 < C 
for all t £ [-00, T 7 ]. 

Step 2: Forward propagation of the regularity 



Using the forward well-posedness of strong solutions stated in Theorem 1.2 (and the discussion that 
follows), we observe that u(t) is a strong solution to the original problem, and so the global attractor 
A is a subset of (H 4 n flg) (O) x flg(fi). 

Step 3: Boundedness of the attractor in (H 4 n ff^) (fi) x H$(Cl) 

In the previous step we have shown that A C (H 4 n Hq) (^) x -^o(^)- However, this does not 
guarantee the boundedness of A in (iJ 4 n Hq) (ft) x Hq(£1). So, using the compactness of the attractor, 
we will follow an additional argument. 

Since for every r £ R, the element u*(r) belongs to a compact set in L 2 (fl) that consists of elements 
from Hq(SI), for every e > there exists a finite set {(f) j} C H^ft) such that we can find indices j± 
and J2 such that 

\\ut(r) - ^jJI + ||ut(r + /i) - 0j 2 || < e. 
Let P(z) be given by (4.9) with the pair w(t) = u(t + h) and u(t) and 

p ji,j 2 ( z ) = - (<f>h ,[ u , v ( z > z )]) ~ (<l>h>[ w > v ( z i z )]) ~ (<t>3i + <f>j 2 , [z,v(u + w,z)]) 
where z(t) = u(t + h) - u(t) = z h (t). It can be easily shown that 

\\P(z)-P n ^(z)\\<tC(R)\\z h (T)\\i (5.4) 

and 

\\Pn, j2 ( z )\\< C{R) (H&Jla + ||^ 2 || 2 )lkV)||L, 

for 77 > 0. So we take 

Bupll^^^ll^CCeJII^MHt, (5.5) 
for 77 > 0. Taking into account (5.4) and (5.5) in (4.8) we see 



sup 

*e[o,r] 



(T(z h ),z' t l ) <C(e,T,R) sup \\z h (r + S )||t„ + e / ||^(r) 

s+t tG[0,T] 



for all s € K with 77 > and T > 0. Now applying the above argument used to prove (5.1) and (5.2), 
we obtain the boundedness of the attractor A in (H 4 n Hq) (SI) x Hq(£1). 

This proves the first part of Theorem |1.4| . 
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5.2 Proof of (ii) in Theorem 1.4 - finite dimensionality 



Using the compactness of the attractor and (5.4)-(5.5), we can write 



T+s 



(.F(*),zt) <C(e)(l + T) sup \\z h {T)\\l_ v + l 

t£[s,s+T] 



s+T 



Ez(t) 



for s € R and 7/ > 0, where z(t) — u(t) — w(t) with (u, Ut) and (w,Wt) from the attractor. Again, 
applying the same procedure from above, we obtain (5.1) for u h = z and s £ I. We then note that 
(5.1) yields 

E z ((m + 1)T) <jE z (mT) + C(A,T)b m , m = 0,1,2, ... 



with < 7 = 7(A) < 1, where 
This yields 



b m = sup \\z(t) 

TE[mT,(m+l)T] 



E z (mT) < 1 E Z (0) + cJ2l m ~ l bi-i 
1=1 



Since 7 < 1, there exist constants Ci,C 2 and a possibly depending on R such that for all t > we 
have 



E z {t) < dE z (0)e- m ( + C 2 sup ||z(r^ 112 



re[Q,t] 



2-J7 



which yields (2.2). Finally, on the strength of Theorem^ applied with B = A, H = D(A 1/2 ) x L 2 (fi), 
"Hi = H 2 ~ v (£l) x {0} we conclude that A has a finite fractal dimension. 
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